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ABSTRACT 

The amount and properties of high-redshift galactic and intergalactic (IG) dust 
are largely unknown, but could be investigated using multi-wavelength photometry 
of high-z objects that have a known intrinsic spectrum. Observations of 7-ray burst 
(GRB) afterglows appear to support the theoretical model of an adiabatic blast wave 
expanding into an external medium. In this model, the synchrotron peak flux is 
independent of frequency, providing a flat spectrum when observed over time, and 
therefore allowing straightforward measurement of the relative attenuation of afterglow 
flux in widely separated bands. Applying this method to dust extinction, we show 
that for a sample of afterglows which have been corrected by galactic extinction, 
comparison between the number counts of peak fluxes in X-ray versus optical can 
provide constraints on an intergalactic component of dust. A similar technique can 
probe the redshift-dependence of extinction in GRB-forming regions without requiring 
an assumed relation between extinction and reddening by the dust. Probing systematic 
changes in extinction with redshift - particularly in IG and/or non-reddening dust - 
is crucial to a proper interpretation of the Type la Supernova Hubble diagram and 
similar observations, and useful in understanding GRB progenitor environments. 

Subject headings: cosmology: observations - dust, extinction - 7-rays: bursts 



Introduction 



There is now substantial evidence that 7-ray bursts (GRBs) originate at cosmological 
distances (e.g. Metzger et al. 1997) from very powerful explosions (e.g. Kulkarni et al. 1998, 
1999). They have been detected at very high redshifts, making them useful probes of the universe 
out to early cosmological epochs. In addition, the 7-ray emission is followed by delayed emission at 
longer wavelengths, from the X-ray to the radio band (Costa et al. 1997; van Paradijs et al. 1997; 
Frail et al. 1997). This afterglow is described reasonably well as synchrotron radiation, emitted 
when a relativistic shell collides with an external medium (Paczyhski & Rhoads 1993; Katz 1994; 
Waxman 1997a, b; Wijers, Rees &: Meszaros 1997; Sari, Piran & Narayan 1998). Afterglows, like 
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their high energy counterpart, are also expected to be detected out to very high redshift. The 
absorption-line systems and the Lya forest visible in their spectra can therefore be used to trace 
the evolution of metallicity in the universe, and to constrain, or possibly measure, the epoch at 
which re-ionization of the universe occurred (e.g. Lamb & Reichart 2000). 

In this Letter we point out other uses that multi-wavelength observations of afterglows might 
have in cosmology. Whereas it is by now clear that GRBs (and consequently their afterglows) are 
far from being standard candles, the afterglows do have a very interesting property: the flux at the 
peak of the synchrotron spectrum is independent of time for adiabatic hydrodynamic evolution 
(Katz 1994; Meszaros k, Rees 1997; Sari, Piran & Narayan 1998), and the adiabatic shock model 
has received robust support from observations (Waxman 1997a,b). Because the frequency of the 
peak flux smoothly decreases with time, this provides a strong theoretical connection between 
observations at various widely-separated frequencies which is lacking in other high-z objects such 
as quasars or galaxies. 

This property makes GRB afterglows well suited to absorption studies of both the GRB 
immediate environment and host galaxy, and of the intervening intergalactic medium (IGM). 
One such possibility is the study of intergalactic dust, which may play a role in observations of 
high-z Supernovae (SNe). In the past few years, observations of Type la SNe by two separate 
groups have revealed a progressive dimming of SNe at high redshift with respect to the predictions 
of a matter-dominated universe or even an open universe with a zero cosmological constant 
(Riess et al. 1998; Perlmutter et al. 1998). This dimming has been interpreted as evidence for 
acceleration in the cosmic expansion, probably caused by a positive cosmological constant. Two 
important systematic effects which must be accounted for in drawing this conclusion are evolution 
in the supernova intrinsic brightness and dust obscuration; both groups have discussed these at 
length. However, their dust corrections, based on reddening, rely on the assumption that dust has 
everywhere the same characteristics as in the Milky Way, which might not be the case. As Aguirre 
(1999) has shown, a scenario in which galaxies expel a significant fraction of their metals in winds 
or by radiation pressure ejection of dust, and in which very small grains are selectively destroyed 
or retained by galaxies, can self-consistently provide a viable candidate for an intergalactic (IG) 
dust component which is more 'grey' (i.e. less reddening) than dust in the Milky Way. Dust 
reddening properties also vary significantly within and between galaxies. 

GRB afterglows are observed in widely-separated bands at which even grey dust will attenuate 
radiation very differently, and, as explained above, a known relation exist between the expected 
fluxes in those bands in the absence of extinction. This Letter shows that, given a sample of 
GRB afterglows for which a reliable correction for the host galaxy absorption has been made, a 
comparison between the number count distribution of peak fluxes in various bands can sensitively 
test the presence of a cosmological distribution of grey dust. Whereas this test is insensitive to 
the choice of model parameters made, such as type of cosmology, shape of the GRB luminosity 
function, redshift evolution of the GRB rate, etc. (because these parameters influence in the same 
way the number counts in all bands, and we are comparing data in different bands), it does rely on 
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the assumption of the constancy of the peak flux in the various bands. In reality, several effects 
can contribute to a departure from the simple, ideal behavior (see e.g. Meszaros, Rees & Wijers 
1998), and we discuss how and to what extent these can be corrected for, so that our test would 
still remain possible. 

Besides probing the existence of IG dust, we show that the same type of multi-wavelength 
study of GRB afterglows can provide a useful probe of the evolution of metallicity (or dust-to-gas 
ratio) in the GRB-forming regions of galaxies. Independent probes of metallicity evolution (i.e. 
not related to GRB sites) can then yield information on the specific environments of GRBs in 
relation to galaxies, helping to elucidate the nature GRB progenitors. Conversely, assuming a 
GRB progenitor type, one could obtain information about dust in the environment of that type of 
object, at very high z. 



2. Number counts of GRB afterglow peak fluxes 

To begin with, we adopt the simplest unbeamed synchrotron model (Waxman 1997a, b). 
Under the assumption that the magnetic field energy density in the shell rest frame is a fraction 
£b of the equipartition value, and that the power-law electrons carry a fraction £ e of the dissipated 
energy, the observed frequency at which the synchrotron spectral intensity of the electrons peaks is 

Mt ) = 2.4 x 10- (i±i)" 2 (il) 2 (£L) 1/2 £ VV 2 Hz , (1) 
where z is the cosmological redshift of the source. For Q = 1 and H = 70 km s" 1 Mpc -1 , the 

1/2 ( £ \ 

observed intensity at v m is F Vm ~ 1 mJy n{ ( or) ^52 for a burst at a redshift z ~ 1. Due 
to the combination of several effects, the afterglow flux falls off with redshift very slowly (Lamb 
& Reichart 2000; Ciardi & Loeb 2000). For an Einstein-De Sitter cosmology, Ciardi & Loeb find 
F Vm cx (1 + z)( p ~ 3 )/ 4 [l — (1 + z) _1 / 2 ]~ 2 , where p is the power-law index of the electron energy 
distribution. Both the GRB and the afterglow observations are well fitted by the value p ~ 2.5 
(Sari, Narayan & Piran 1998; Kumar & Piran 1999). We fix the constant of proportionality so 
that F Um = 1 mJy at z = 1. In absence of absorption, the peak luminosity of a given burst is then 
given in terms of the luminosity distance Dl(z) by L Vm = \^D\jyZ)F Vm j (\ + z) (see e.g. Hogg 
1999 and references therein). At each redshift, we allow for a scatter in the burst luminosity by 
assuming a probability distribution that is log-normal in L Um . The mean of this distribution is 
chosen so that {F Um (z)) evolves according to Ciardi & Loeb's relation. We choose a = 0.5 in our 
computations. However, as explained in the following, our results are not very sensitive to the 
particular choice of parameters. 

The number (per unit time) of GRB afterglows with peak flux above a given threshold S is 
given by 

N(> S) = / dzR(z) -f \ dL Um P(L Um ,z) . (2) 

Jo dz Jl^U^s) 
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Here R(z) is the GRB formation rate, which we assume to be proportional to the star formation rate 
(SFR), as given by the most recent observations of Steidel et al. (1999), and z max is the redshift limit 
of the sample. The comoving volume is given by dV c = 4(c/# ) 3 (l + z- ^/T+^) 2 (l + z)- 7 l 2 dndz 
(e.g. Hogg 1999) in our adopted cosmology. The minimum luminosity that a burst at redshift z 
must have in order to be detected with flux > S is given by 



Here r a b s is the optical depth to photoelectric absorption, for which there is a contribution from 
the host galaxy itself and a contribution from the intergalactic medium (IGM). X-rays are not 
affected by photoabsorption, if observed at high enough frequency (see below). Optical light is 
unabsorbed if observed below the Lya resonance frequency v a {z) = 2.47 x 10 15 /(1 + z) Hz, so the 
y-band {y = 5.4 x 10 14 Hz) and longer wavelengths will not be affected by photoabsorption if one 
only considers bursts with redshifts up to z ~ 3.5. We will show our results by assuming a sample 
with z max =3. 

The optical depth of dust at redshift z is indicated by T^ ust in Equation (g). We assume 
T dust( z ) = T gai( z ) + T grcy(z) , and consider several possibilities. 

For the grey dust component, we take the model as in Aguirre (1999): dust is composed 
of equal masses of graphite and silicate spherical grains of radii a, with a grain-size distribution 
dN(a)/da oc a~ 3 ' 5 , 0.1 /um < a < 0.25 ^m. This dust has a very fiat extinction curve blueward of the 
V band. We take a dust density fidust(£ = 0) = 4.5 x 10~ 5 , and assume that rid U st(-2) oc / dtSFR, 
where the SFR is the same as that assumed for the GRB rate. Using publicly available extinction 
data (see Laor & Draine 1993) these assumptions give a cosmologically 'interesting' extinction of 
~ 0.2 mag to z = 0.5 in the observed B-band. 

For the host galaxy absorption we assume the extinction law to be the same as in our Galaxy, 
and model, for an observed wavelength A, 



with £(A) fitted as in Pei (1992). We then consider two types of evolutionary scenarios for the 
dust-to-gas ratio k(z); model (1): k(z) = const = k(0), with k(0) = 0.78 as in the Galaxy; model 
(2): k(z) traces the integrated star- formation rate. These should represent the most extreme 
possibly evolutionary paths for the dust-to-gas ratio: the second effectively assumes instantaneous 
mixing of metals throughout all cosmic gas, whereas the first effectively assumes either that cosmic 
metallicity does not evolve, or alternatively that metals are so poorly mixed that while cosmic 
metallicity evolves, the metallicity near the GRBs does not. Model (1) could better represent 
a scenario in which GRBs are associated with coalescent compact objects (e.g. Eichler et al. 
1989, Narayan, Paczyhski & Piran 1992), and therefore they would not be expected to occur in 
particularly dense or dusty regions. On the other hand, if GRBs are associated with the death of 
massive stars (e.g. Woosley 1993, Paczyhski 1998), then they are likely to occur in star- forming 
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(1 + z) exp{-[r dust (z) + r abs (^)]} ' 
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regions. The mean metallicity expected for such regions is substantially higher than the cosmic 
mean. This is a situation which could be more closely described by our dust model (2) 

In Equation (H), JVh(z) is the average column density intercepted by the GRB photons as 
they escape a host galaxy at a redshift z. To estimate this function, we have performed a Monte 
Carlo calculation using the model of Ciardi & Loeb (2000) to obtain the probability distribution 
Nji(M,z) for GRBs occurring (with spatial distribution proportional to the gas density squared) 
in the disks of spirals embedded in halos of mass M, combined with the Press-Schechter formalism 
for the distribution of dN(M, z)/dM. The resulting distributions are normal in log(Nff). As z 
runs from to 4, the distribution's center increases roughly linearly from 21.27 to 21.48, and its 
width decreases linearly from 0.915 to 0.73.0 This should be regarded not as an accurate and 
detailed model of the expected column densities (for example it does not include the evolution of 
the galaxies' gas fraction), but rather as a rough check on the possible systematic increase of GRB 
host-galaxy extinction. Our two models, based on this estimate, bracket the more detailed models 
of extinction evolution of Calzetti & Heckman (1999) and Totani & Kobayashi (1999) which apply 
to the normal stellar population in galaxy disks. 

3. The effects of dust on GRB number counts 

We present results for number count distributions with an observed peak flux in X-ray, V, 
and J bands. Observations in the hard X-ray band are necessary for a proper comparison which is 
not biased by photoelectric absorption within the host galaxy. This should not be a problem for a 
5 KeV band at which Swift can observe: Ghisellini et al. (1999) find that, at this energy, a burst 
at z = in a medium with solar metallicity will be unaffected by absorption for column densities 
Nn <^ 10 23 cm~ 2 . For a burst at z ~ 2 in a solar metallicity environment, the constraint becomes 
Xh <^ 10 24 cm -2 . Evidence for an unusually high column density to a burst can be found by the 
time-dependence of the UV flux (Perna & Loeb 1998), K a fluorescent line emission (Ghisellini et 
al. 1999), and, clearly, heavy reddening of the optical emission. For a proper analysis such bursts 
should not be included in the sample. 

To probe the presence and type of dust extinction, the observed optical-near IR bands seem 
the most appropriate. As explained in §2, the V band and longer wavelengths are unaffected by 
photoelectric absorption in the IGM up to z ~ 3.5. Therefore, if the sample is limited to redshifts 
lower than this, the on-board Swift V photometry could be used to determine the peak V flux 
(although ground-based NIR data may be necessary for obtaining the rest-frame UBV photometry 



2 Note that, if GRBs are associated with star-formation regions, then also the density in their surrounding medium 
is expected to be on average higher. However, extinction is only sensitive to the product of column density and 
metallicity, and therefore we parameterize our models with only one of them, that is the dust-to-gas ratio. 

3 For (GRB density) oc (gas density), the mean values are similar but about 0.25 dex smaller. 
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necessary for accurate de-reddening). 

Figure 1 shows a comparison between number counts in X-ray, V, and J bands for a sample 
of GRBs for which perfect corrections for Galactic and GRB host-galaxy extinction have been 
made, while the grey component (which would not be accounted for by UV/optical de-reddening) 
has not been subtracted out. Figure 2 shows the same curves, where host galaxy extinction - 
quantified by models (1) and (2) of §2 - has not been corrected. In these figures, the horizontal 
offset of two curves indicates the amount of extinction, whereas the vertical change in that offset 
shows the evolution of extinction with increasing redshift. 

In a realistic case, one could hope to correct for host galaxy extinction only imperfectly: 
there will be errors in the photometry and .fT-correction of the rest-frame UVB fluxes, errors 
in determining the appropriate extinction law (Ry, see Cardelli, Clayton & Mathis 1989), and 
variations in the extinction law even for a given Ry. The extinction may also be time-dependent 
(Waxman & Draine, 1999). However, we find that even if the correction leaves a 'residual', the 
statistical difference is still readily measurable. We performed a Monte Carlo simulation of data 
drawn from the X-ray distribution and of data drawn from the T^-band distribution corrected for 
galactic extinction but not for a grey intergalactic component (as in Figure 1). To each flux was 
assigned a fractional error that was randomly generated from a Gaussian distribution of width 
(Terr- We took <7 err = 0.5, which corresponds to about half to 1 magnitude error in extinction 
correction, larger than the quoted errors on A v for GRBs observed so far (see e.g. Bloom et al. 
1998, Vreeswik et al. 1999). By performing a Kolmogorov-Smirnov test (see e.g. Press et al. 1992) 
with the two sets of simulated data, we found that with N ^ 70 bursts the two distributions can 
be distinguished with a confidence level a p> 3. 

Note that, if the host galaxy extinction is decreasing with z (as in our model (1), or as in 
Calzetti & Heckman 1999 and Totani & Kobayashi 1999), then galactic extinction would never 
mimic the effects of IG grey dust. On the other hand, if the likely scenario is as in model (2), then 
our ability to distinguish IG dust from increasing host galaxy extinction would depend on the 
accuracy with which host galaxy dust correction can be made. 

A very clear distinction can also be made between the distributions described by model (1) 
and (2) in Figure 2. A random simulation as described above showed that with N ^ 20 bursts the 
two y-band distributions in the case where there is no grey dust component can be distinguished 
with a confidence level a ^ 3. Similar results hold for the corresponding distributions which also 
include grey dust. In these simulations we have assumed the same <T er r as above. However, when 
no correction for galactic extinction is required, this error is likely to be much smaller. This study 
is therefore very sensitive to the conditions (i.e. metallicity and density) in the environment of 
GRBs. Combining such analysis with independent probes of the metallicity and density evolution 
of the medium in galaxies and in the IGMj], significant information can therefore be obtained 



4 Spectra of the GRBs themselves could also be used to such purposes, as shown by Lamb & Reichart (2000). 
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about the location of GRBs within galaxies. This is particularly relevant in order to distinguish 
scenarios in which GRBs are due to coalescence of two compact objects with respect to scenarios 
in which GRBs are associated with the death of massive stars. Since reddening and extinction 
are measured rather independently, the technique can also give information about the reddening 
properties of dust in GRB progenitor environments. 

We need to emphasize once more that our method is rather insensitive to the choice of 
model parameters made, such as type of cosmology, shape of the GRB luminosity function, 
GRB evolution with redshift, etc. This is because these parameters influence in the same way 
the number counts in all bands, and our method is based on a comparison of data in different 
bands. What can change (due for example to a different distribution of the bursts with the 
redshift) are statistical quantities such as the number of bursts required to distinguish between 
two models with a given confidence level, but we found the difference in this number required 
for various models to be rather marginal. The main uncertainty of this method lies in the fact, 
already discussed in §1, that the perfect, adiabatic model with a constant peak flux is a clear 
approximation that does not take into account several effects which might cause a departure from 
this simple behavior. Meszaros et al. found that a number of factors (such as angular anisotropy 
of the fireball, properties of the environment, later re-energization of the afterglow, etc.) can lead 
to either an increase or to a decrease of the peak flux with time, depending on what the particular 
conditions (such as for example density gradients in the external medium) are for a given burst. 
This means that, most likely, such disturbances should not give rise to systematic effects, but 
rather would lead to variations equivalent to having a larger scatter in the errors on the fluxes. 
If this is indeed the case, the effects produced by grey dust should not be easily washed out: a 
random simulation of data with errors in peak fluxes drawn from a distribution with a crr =2 (i.e. a 
fractional error of 200%) showed that the distributions in Figure 1 can still be distinguished with 
a confidence level a ^ 3 if the number of afterglows is ^ 200. Clearly, the situation would become 
more difficult if these disturbing effects, instead of being stochastic, had a trend which would 
systematically reduce the emission in one (or more) bands with respect to the others. However, 
unless these hypothetical systematics have also an evolution with redshift, it should be possible to 
correct for them by measuring ratios between peak fluxes in various bands for a sample of bursts 
at low z (which would only be marginally affected by IG dust). A systematic decrease in the flux 
at peak frequency could also mimic the difference between a dusty and a relatively dust-free GRB 
environment, but again only at one redshift; evolution in the GRB environment should still be 
measurable. Such a study, to be statistically significant, would need a large number of observed 
afterglows. This should be possible with Swift, which will observe ~ 300 bursts per year, and 
will monitor their emission from the X-ray to the V band (and other instruments will very likely 
obtain NIR photometry of the afterglows). It will be an ideal tool for such a study. 
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4. Conclusions 

One of the outstanding problems in observational cosmology is to understand the type and 
amount of dust existing at high-z, whether outside galaxies or within them. Traditional methods 
of estimating extinction require a known dust reddening curve and an accurate model for the 
intrinsic UV/optical spectrum of the observed object, both of which may be rather uncertain; very 
few existing techniques allow the measurement of absolute extinction. 

In this Letter, we have proposed GRB afterglows as probes of dust in cosmology. Like 
QSOs and galaxies they can be observed to very high z. In addition, there is a well-based 
theoretical spectrum which is flat between X-ray and near infra-red frequencies. This allows very 
straightforward estimates of dust extinction independent of its reddening properties (though the 
reddening would be observable and useful) in GRB environments or in the intergalactic medium; 
the former can help help determine the properties of GRB progenitor environments (including the 
dust properties), while the latter is crucial for careful interpretation of the Type la Supernova 
Hubble diagram. 

While relying on the assumption that afterglow evolution is well described by the simplest 
adiabatic blast model, with a large number of afterglows (such as Swift will provide) , the technique 
should be applicable even if deviations from the theory occur. 
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Fig. 1. — Cumulative number count distribution of GRB afterglow peak fluxes in the X, J, and V 
bands for a sample corrected by galactic extinction but not by the extinction due to an intergalactic 
grey dust component. 
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Fig. 2. — Cumulative number count distribution of GRB afterglow peak fluxes in the X, and V 
bands (upper panel) and the X and J band (lower panel). In both panels, the dust models (1) and 
(2) refer to the host galaxy dust models described in the text. For each model, the upper curve has 
only the component from galactic extinction, while the lower curve includes also an intergalactic 
grey dust component. 



